Weak superconducting pairing and a single isotropic energy gap in stoichiometric LiFeAs 
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We report superconducting (SC) properties of stoichiometric LiFeAs (T c = 17K) studied by small-angle 
neutron scattering (SANS) and angle-resolved photoemission (ARPES) . Although the vortex lattice exhibits 
no long-range order, well-defined SANS rocking curves indicate better ordering than in chemically doped 
122-compounds. The London penetration depth A ab (0) = 210 ± 20 nm, determined from the magnetic 
field dependence of the form factor, is compared to that calculated from the ARPES band structure with no 
adjustable parameters. The temperature dependence of X ah is best described by a single isotropic SC gap 
A = 3.0 ± 0.2 mey which agrees with the ARPES value of a£ rpes = 3.1 ± 0.3 meV and corresponds to the 
ratio 2A//c B T c = 4.1 ± 0.3, approaching the weak-coupling limit predicted by the BCS theory. This classifies 
LiFeAs as a weakly coupled single-gap superconductor, similar to conventional metals. 

PACS numbers: 74.70.Xa 61.05.fg 74.25.-q 74.25.Uv 



In many of the recently discovered Fe-based superconduc- 
tors (SC) [1, 2], a transition to the SC state is induced by 
chemical doping of a parent compound that at ambient con- 
ditions does not exhibit SC in its stoichiometric composition 
even at the lowest temperatures. Among the few known ex- 
ceptions, the present record holder for the SC transition tem- 
perature, T c , is the stoichiometric LiFeAs (T c < 18 K) [3-5]. 
Others are low-T c superconductors NaFeAs (T c = 9K) [6], 
FeSe (T c = 9K) [7], LaFePO (T c = 6K) [8, 9], and KFe 2 As 2 
(T c = 3.8 K) [10]. The electronic structure of LiFeAs is quasi 
two-dimensional (2D) [11] and supports superconductivity 
in the absence of any notable Fermi surface (FS) nesting or 
static magnetism [12]. However, the presence of normal- 
state antiferromagnetic fluctuations has been suggested by 
75 As NMR measurements [13]. Together with the weakness 
of the electron-phonon coupling predicted by the density 
functional theory [14], this suggests that the SC pairing in 
this structurally simple compound possibly has the same 
magnetic origin as in higher- T c iron pnictides [15-17]. On 
the other hand, arguments advocating the phonon mech- 
anism have also been raised recently [18]. Therefore, to 
pinpoint the SC mechanism with certainty, details of the SC 
pairing symmetry and the coupling strength are required. 

In a number of recent studies [19-22], it was shown that 
doped iron arsenide superconductors are characterized by 
strong pinning of magnetic flux lines that precludes the 
formation of an ordered Abrikosov lattice. The role of the 
pinning centers can be played by magnetic/structural do- 
mains in the underdoped samples [21], by the dopant atoms 
themselves, such as Co or Ni, at higher doping levels [22], or 
by the electronic inhomogeneities that result from phase sep- 
aration in some hole-doped 122-systems [23]. This served 
as our motivation to study the magnetic field penetration 



in a single crystal of stoichiometric LiFeAs, which possesses 
a non-magnetic ground state with tetragonal crystal sym- 
metry, thus excluding all of the above-mentioned strong 
pinning mechanisms from consideration. In the following, 
we will compare these results with ARPES measurements of 
the electronic structure to establish the microscopic origin 
of the measured quantities. 

For our SANS experiments, we used a large single crys- 
tal with a volume of~10xl0x 0.4 mm 3 . It was grown 
by the self-flux method and characterized as described in 
Ref. 12. The value of T c , measured on a smaller piece of 
the same sample, was ~ 1 7 K, as shown in Fig. 1 (a) . Mag- 
netization data are corrected for demagnetization by an 
ellipsoid approximation [24] . In order to avoid any expo- 
sure of the sample to air during mounting into the cryostat, 
it was placed inside a single -crystalline Si box [Fig. 1(b)]. 
The box was sealed with a thin rectangular lid prepared 
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Fig. 1 (color online) . (a) Magnetic susceptibility of LiFeAs, mea- 
sured upon warming after cooling in magnetic field (FC) and 
in zero field (ZFC) . (b) Photo of the sample prepared for SANS 
measurements inside the single-crystalline silicon box (see text) . 
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Fig. 2 (color online). (a, b) SANS diffraction patterns measured 
at H = 0.25 and 0.5 T, respectively. The 0.25 T data are summed 
up over the rocking curve from —8° to 8°, whereas the 0.5 T data 
are shown for the zero rocking angle. Both datasets are smoothed 
with a 3-pixel FWHM Gaussian filter, (c) Angle-averaged diffracted 
intensity as a function of momentum transfer |q|, measured at 
different magnetic fields between 0.25 T and 1.5 T. For clarity the 
zero line of each curve is offset from the one below it. Vertical 
arrows show the expected peak positions for a perfect triangular 
lattice q A . Solid lines are Gaussian fits, (d) Averaged intensities 
on the left (❖) and right (o) sides of the ring as functions of the 
rocking angle, measured at H — 0.25 T. 



from a Si wafer that was glued on top with a small amount 
of low-temperature varnish inside the glove box with Ar 
atmosphere. Because of the low absorption and incoherent 
neutron-scattering cross-sections of Si, such box is essen- 
tially transparent to neutrons and has a negligible effect 
on the background. To mark the position of the sample 
inside the sealed box with respect to the neutron beam, two 
stripes of highly neutron-absorbent Cd were placed above 
and below it. 

The SANS experiment was carried out using the Dll 
instrument at the Institut Laue-Langevin (ILL) in Grenoble. 
We used the usual experimental geometry, with the applied 
field approximately parallel to the incident cold-neutron 
beam, characterized by the wavelength A = 8 A and the 
full width at half maximum (FWHM) wavelength spread 
A A/ A = 10%. The diffracted neutrons were collected by 
a 2D multidetector placed behind the sample. The vortex 
lattice (VL) was prepared in the sample by applying the 
desired field above T c , and subsequent field-cooling to 2 K. 
In all cases, background measurements were carried out 
at 20 K (above T c ) and subtracted from the field-cooled 
foreground measurements. 

Two representative diffraction patterns measured in mag- 



netic fields H = 0.25 and 0.5 T are shown in Fig. 2 (a, b). 
A distinct ring of scattering is seen in both panels. The ab- 
sence of distinct Bragg peaks indicates a lack of long-range 
orientational order of the VL. The 0.25 T image is a sum of 
10 individual diffraction patterns measured at different rock- 
ing angles between —8° to 8°, whereas the 0.5 T image was 
measured at zero rocking angle only. At higher fields, the 
scattering intensity is decreased, but a peak can be seen in 
the angle-averaged data up to H = 1.5 T, as shown in panel 
(c) . For all fields, the fitted peak position agrees within the 
experimental error with the expected radius of the diffrac- 
tion ring q A = 4n/a A ^, calculated for a perfect triangular 

VL with lattice parameter a A = yj 2<fi /HV3, where O is 
the magnetic flux quantum. Calculated positions of q A are 
shown in the figure by small vertical arrows. 

So far, these results are similar to those reported for 
electron-doped 122-compounds [19, 22]. The first essen- 
tial difference is illustrated in Fig. 2 (d), which shows rock- 
ing curves with half-width at half-maximum (HWHM) of 
3.2 =b 0.5°, measured at 0.25 T. They represent the angular 
dependence of the diffracted intensity on the left and right 
sides of the diffraction ring. Despite the clear intensity vari- 
ation, disorder-induced pinning causes the rocking curves 
to remain broad, and using the HWHM we estimate the lon- 
gitudinal correlation length of the vortices as f y w 0.8 /im. 

Conversely, in previous SANS experiments on both un- 
derdoped [19] and overdoped [22] BaFe 2 _ x Co x As 2 , the 
rocking curves were much broader and extended beyond 
the measurable range. Our observations therefore indi- 
cate an improvement in the longitudinal VL ordering and 
a decrease of the typical pinning forces in the absence of 
chemical dopants. But the fact that we did not observe any 
long-range orientational VL order even after oscillating the 
field value by 1 and 10% during in-field cooling suggests 
that the pinning in our sample is still not negligible. 

Now let us turn to the quantitative determination of some 
important SC properties of LiFeAs. The integrated intensity 
J corresponding to 1 /6 of the diffraction ring (one Bragg 
spot of a triangular VL), obtained from the rocking curve, is 
proportional to the modulus squared of the VL form factor 
F(q,T) [25], i.e. the Fourier transform of the 2D magnetic- 
flux modulation within the sample. As our measured inten- 
sity mainly originates from the first-order Bragg spots at a 
distance q&q A from the origin, we limit our considerations 
to the first order VL form factor, 

/ = 27ry<D( r /4) 2 A 2 0- 2 q- 1 |F(q,r)| 2 . (1) 

Here V is the sample volume, 3> is the neutron flux density, 
and y is the magnetic moment of the neutron in nuclear 
magnetons. By varying the magnetic field and, consequently, 
q A , we can thus study the q -dependence of the form factor, 
assuming that the rocking curve width is both field- and 
temperature-independent. 

In the simplest Ginzburg-Landau model valid for super- 
conductors with large k = A/£ » 1, and in small magnetic 
fields H <C H c2 , the form factor can be expressed in terms of 
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Fig. 3. (a) Magnetic field dependence of the VL form factor at 
T = 2 K, fitted to Eq. (2). (b) Temperature dependence of |F(q,T)| 2 
at H = 0.25 T. The vertical axis is scaled to the value of A -4 (0) 
that resulted from the form-factor fit in panel (a). For comparison, 
A _4 (T) for the d-wave gap is shown by the dashed line. 



the temperature-dependent penetration depth X ah {T) and 
the SC coherence length £ a b(T) [26] (for brevity, the index 
ab will be skipped), 

Here K^g) is the modified Bessel function of the second 
kind [38] . Because the band structure of LiFeAs can poten- 
tially support two-gap SC [12], we resort to the following 
analytical two-gap model to describe the temperature de- 
pendence of the penetration depth, 
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where the function M is an accurate approximation [27] for 
the temperature-dependent quasiparticle effects, and the 
constants I 12 depend only on the band structure and can 
be expressed as simple integrals over the FS [27, 28]. The 
T-dependence of the SC gap is approximated by [29] 



A(T) = A tanh^y/ TjT- 



(4) 



Consequently, because in the weak-coupling limit the SC 
coherence length is inversely proportional to the energy gap 
[30], the temperature evolution of E, will be given by 



tarihC^y/TjT-l 



oc 1/A(T). (5) 



Substituting expressions (3) through (5) into Eq. (2), we 
obtain the final formula for fitting the SANS data. 

Following the methodology of Ref. 31, we start with the 
magnetic field dependence of the low-temperature form 
factor that is shown in Fig. 3 (a). Because at T = 2K the 
values of £ and A can be considered equal to their zero- 
temperature limits, Eq. (2) can be applied directly to the 
field-dependent data (solid line in the figure) to extract the 
values of A ab (0) = 210 ± 20 nm and £ ab (0) = 7 ± 2nm 
(k = 29 ± 7). The obtained value of X ah agrees with the 



results of a muon-spin rotation (/iSR) measurement [32], 
which yielded A ab (0) = 195 and 244 nm for two samples of 
Li 1+5 FeAs with T c = 16 and 12 K, respectively. Our value of 
£ a b(0), however, is likely to be overestimated with respect to 
that obtained from upper-critical-field measurements [33], 
which is £ab(0) = <J 4> /2nH^ 2 & 2 — 4nm. Such overesti- 
mation can result either from a finite-width T c distribution 
in our large sample, or from the field-induced disorder of 
the VL expected in the Bragg glass model [34] . 

Next, we turn to the temperature evolution of the form 
factor. As follows from Eq. (2), for A(T)q » 1, F(q,T) oc 
A~ 2 (T). The scattered intensity therefore scales oc A -4 . In 
Fig. 3 (b), the measured integrated intensity is plotted as a 
function of temperature, and the vertical axis is scaled to 
the value of 1/A 4 (0) that resulted from the low-temperature 
fit of the form factor. By fixing A(0) and £(0) to the values 
found previously, we can now fit the two SC gaps A 1 2 (0) 
and the coefficients using |F(q,T)| 2 as the fitting func- 
tion. It turns out that independently of the parameter ini- 
tialization, the fit converges to a single value of the gap 
A 1 = A 2 = 3.0 ± 0.2 meV This value of the energy gap 
corresponds to the ratio 2A/k B T c = 4.1 ± 0.3, approaching 
the weak-coupling limit of 3.53 predicted by the BCS theory 
of conventional superconductivity [35]. For comparison, 
A~ 4 (T) corresponding to a d-wave gap is also shown in the 
same figure, producing a poor fit. This essentially excludes 
the possibility of two-gap SC or gap nodes in LiFeAs. 

Now we compare these results with those of ARPES, to 
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Fig. 4 (color online) . ARPES spectra of LiFeAs measured on the 
double-walled electron-like M-barrel in the SC (a) and normal (b) 
states, (c) The integrated energy distribution curves (IEDCs) of the 
same spectra, (d) The low-temperature IEDC after normalization, 
fitted to the Dynes function. 
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establish their relationship to the microscopic electronic 
properties, such as band dispersion and the SC gap. An anal- 
ysis of the leading edge shift along the FS contours implies 
an isotropic gap for every FS sheet [12] . To quantify the low- 
temperature gap value A , we employed the Dynes function 
fitting procedure [36] to the ARPES spectra measured on 
the double-walled electron-like M-barrel [Fig. 4 (a, b)]. The 
energy distribution curves integrated in a wide momentum 
window along the FS radius (IEDCs), measured in the SC 
state below 1 K and in the normal state at 23 K, are shown in 
Fig. 4 (c) . In order to reveal the true shape of the spectrum 
in the SC state, the low-temperature IEDC was normalized 
by the Fermi-function-corrected normal state spectrum, as 
shown in Fig. 4 (d). The good quality of the Dynes -function 
fit confirms the robustness of such normalization. The re- 
sulting low-temperature value of Aq RPES = 3.1 ± 0.3 meV 
is in perfect agreement with that extracted above from the 
temperature dependence of X ah . 

The knowledge of the band dispersion together with the 
SC gap allows the calculation of macroscopic properties in 
the SC state with no adjustable parameters. The superfluid 
density at T — > is proportional to the integral of Fermi 
velocity v F along the FS perimeter [27, 28, 37], and in the 
clean limit, 
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where e , h, e, c are physical constants, and L c is the c-axis 
lattice parameter. Although the FS of LiFeAs consists of sev- 
eral electron- and hole-like sheets [12], for the evaluation 
of the integral (6) the renormalized Fermi velocity, extracted 
from ARPES data, can be well approximated by its average 
value of h{v ¥ ) = 0.3 ± 0.03 eVA. For the experimental LiFeAs 
band structure, this formula yields A^ PES = 172 ± 20 nm, 
which is only slightly lower than our directly measured 
value. Similarly, at T — » 0, the BCS coherence length is 
proportional to the ratio of Fermi velocity to gap magnitude, 
^arpes^q) = fr( VF )/ u A [35], which equals 3.2 ± 0.4 nm in 
our case. This corresponds to the upper critical field H^ 2 = 
32 ± 8 T, in agreement with direct measurements [33] . 

In summary, we have evaluated several important SC pa- 
rameters of LiFeAs from two complementary experiments. 
We have demonstrated that its order parameter is isotropic 
and in contrast to the higher- T c ferropnictides [27] is char- 
acterized by a single SC gap A = 3.0 ± 0.2 meV This value 
is close to the BCS limit of 1.76 k B T c , which indicates that 
LiFeAs is a weakly coupled single-gap superconductor, simi- 
lar to conventional metals. 
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